Propriedades térmicas 
dos gases 1: equações 


de estado 


META 

Generalizar o formalismo de mecânica 
estatística clássica visto na aula passada 
para obtenção das equações de estado de 


gases reais. 


OBJETIVOS 

Ao fim da aula os alunos deverão ser 
capazes de: 

Entender a modelagem de gases reais. 
Compreender o conceito da função de 
grande partição. 

Obter equações de estado usando esta 
função. 

Resolver problemas envolvendo estes 


conceitos. 


PRÉ-REQUISITOS 
Aulas anteriores e cálculo diferencial e 


integral. 
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6.1 Introdução 


Prezado aluno, apresentamos na aula anterior conceitos de me- 
cânica estatística clássica. Que tal agora usarmos estes conceitos 
para estudarmos sistemas clássicos? É justamente isto que iremos 
fazer nas próximas duas aulas, explorando propriedades térmicas 
(termodinâmica) dos gases. 

Nesta aula iremos estudar as equações de estado (rever sec. 


2.2) para o gás ideal e para gases reais (aproximações). 


6.2 Gás ideal 


Sabemos que a entropia de equilíbrio estatístico é (rever ativ. 
5.1) 
U 


Z 
S=m+hksNlIn+hkoN. 


Considerando o número de partículas constante, dN = 0, e difer- 


enciando esta equação, temos 


dU UU 


Agora podemos multiplicar esta equação por T e lembrar da re- 


lação entre calor e entropia dQ = TdS. Com isso, obtemos 
U 
dQ = dU — qir + NkpTdln Z. 


Pela primeira lei da termodinâmica temos dQ = dU + PdV. Logo, 
encontramos 


NksTdin Z = SdT + Pav. 


Portanto, concluímos que 


U=NkT E mZ, (6.1) 
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ô 
P=NkpT sz (6.2) 


Note que a eq. (6.1) é equivalente à eq. (5.23). 
Na aula passada, calculamos a função de partição de um gás 


ideal [eg. (5.27)], que é 
Z = V(2amksT)?, 
e mostramos que 
3 3 
U = NkpT, ou U=SNRT. (6.3) 
Usando a eq. (6.2), encontramos 
PV = NkpT, ou PV=NRT, (64) 


que é justamente a lei do gás ideal. As eq. (6.3) e (6.4) são as 


equações de estado para um gás ideal. 


6.3 Gases reais 


No tratamento de gases reais clássicos, devemos considerar as 
interações intermoleculares e as dimensões finitas das moléculas. 
Como as forças entre as moléculas são de alcance curto, quanto 
mais diluído for o gás real, mas próximas serão suas propriedades 
das de um gás ideal. A concentração (inverso da diluição) do gás 
pode ser quantificada pelo número de partículas por volume N/V. 
Sendo assim, podemos considerar a seguinte aproximação para um 
gás real muito diluído 

NkBT 


Px S 
V 


Podemos imaginar que quanto menor a diluição (maior a concen- 


tração), temos que considerar termos de ordem superior de N/V, 


“A denominação virial 
é referente ao teo- 
rema do virial visto 
em mecânica de um 
sistema de muitas 
partículas interagentes 


[E]. 


—s, 
E ( 


(108) 
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ou seja, devemos levar em conta a seguinte expansão polinomial 


NkBT N? Nº Nº 
SE a O ua (6.5) 


P 
V V2 Vs 


sendo esta a equação de estado de um gásreal. 4, B,e €C dependem 
de cada gás e são chamados, respectivamente, de segundo, terceiro, 


quarto, ..., coeficientes do virial*. 


Até aqui, fizemos uma abordagem qualitativa. No entanto, 
usando os métodos de mecânica estatística, é possível estimar os 


coeficiente do virial para alguns casos particulares. 


6.3.1 Energia total de um sistema de partículas iso- 
lado 


Em geral, a energia total de um sistema isolado de partículas in- 


teragentes é 


Simplificando a notação, temos 


2 
Ki =K(pi) = + 


“9m 


sendo a energia cinética da partícula 1 com momento p; e 


sendo a energia potencial ocasionada pela interação entre a partícula 
i, com posição 7;, e a partícula j, com posição 7;. Com isso, es- 


crevemos 


E = SM E SU: (6.6) 


i<j 
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O somatório em à < j é para evitar termos de autointeração (à £ j) 
e a contagem de termos repetidos, pois a interação da partícula à 


com a j é a mesma da j com a i, ou seja, U;; = Uj;. Sendo assim, 


> ui; = Ur + Uis 
i<j 


S 


Perceba que a energia total do sistema não é a soma da ener- 
gia de cada partícula isoladamente, pois a energia potencial é um 
termo que depende das partículas dois a dois. Sendo assim, temos 
que adaptar o formalismo estatístico visto na aula passada para 


levar em conta esta situação. 


6.3.2 Função de grande partição 


Como visto na aula passada, para um sistema de partículas 
idênticas e não interagentes, substituindo a eq. (5.18) na eq. 


(5.17), obtemos 
Z= / de / dir / dêpo[e — e(r,p)| exp(— Be). 
0 


Integrando na energia e usando as propriedades da ô de Dirac, 


temos 
Zi : dr / Epexpl-Be(?, 7], (6.7) 


ou seja, a função de partição da partícula é a integral em todo 
o espaço de fase do peso de Boltzmann, exp[-Be(r,D)], e e(7,D) 
é a energia da partícula. No entanto, agora precisamos tratar o 
sistema como um todo, ou seja, não podemos tratar suas partículas 
isoladamente. Por isso, devemos também generalizar a eq. (6.7). 


Para isso, introduzimos a função de grande partição, definida por 


1 E E s, 5 Em 
Z= m/ en [amo [dry | Eixel-s2(a TE) 
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Perceba que agora integramos o espaço de fase de todo o sistema 
(todas as partículas) e por isso, a energia no peso de Boltzmann 
deve ser a energia total do sistema. A divisão por N! continua 
sendo a correção de Gibbs, a qual é necessária para que a conexão 
entre mecânica estatística e termodinâmica seja coerente. Pode- 


mos re-escrever a eg. acima com uma notação mais compacta 


N 
1 =; > 
= ai [am [ato exp(—5E), (6.8) 
Ddr 


onde fica implícito que E é a energia total do sistema, a qual 
depende dos momentos e posições de todas as partículas. 

Perceba que se as partículas não interagem, a energia total 
do sistema é a soma da energia de cada partícula, ou seja, E = 


ei+::-+en. Neste caso, temos 
1 N N 
Dn Ni [er [as exp [-254) 
j=1 i=1 
1 N N 
7 [é [as [exp(-se;) 
j=1 i=1 
1 N 
= mi [dm [ dojem(-s6) 
j=1 


f N 
= ml z (6.9) 
j=1 


onde Z; é a função de partição da j-ésima partícula. Se consider- 
armos que as partículas são idênticas, devemos fazer Z; = Z, para 


todo 3. Sendo assim, obtemos 
A qe (6.10) 


Baseando-se na eq (6.10), devemos reformular as egs. (6.1) e 


(1104 
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(6.2) da seguinte forma 


2 

2 

U=knT plz, (6.11) 

PER nE (6.12) 
Gui Té 


No caso geral de partículas interagentes, podemos substituir a 


eq. (6.6) na eq. (6.8) e, com isso, escrevemos 


2 = a II / é fes exp | — APR + Ui; , 


i<j 


E nm 1/0 exp (= 8Kk) (1/4 a) sata 
ad Ei (1/4 a) [exp (au) 


i<j 


= der à (ya / 85) Mest — Bis) 


1<] 


onde Zgy é a função de partição do gás ideal, pois cada integral 
em d?p; corresponde à (2xmksBT)*/2 = Ze1/V. Como todas as N 
integrais, possuem o mesmo valor, o produtório em k faz com que 
apareça o expoente N. Agora, podemos recorrer à eq. (6.10) e 


escrever 


Z = ZgZim (6.13) 


onde Zey é a função de grande partição do gás ideal, que é 


(CamksT)N'2, (6.14) 


Zi = [7x1 / 4) [exp (su). (6.15) 


(113, 
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Figura 6.1: Comportamento da energia potencial intermolecular em função 


da distância das moléculas. (Figura adaptada da ref. [2].) 


Para conhecermos as propriedades termodinâmica de um gás ideal, 
precisamos encontrar sua função de grande partição. Para isso, de- 
vemos conhecer a expressão da energia potencial entre suas partícu- 


las U;;, pois assim, seremos capazes de desenvolver Zint. 


6.3.3 Segundo coeficiente do virial 


Considere a expansão do peso de Boltzmann em série de potên- 


cias 
exp (—BUij) = 1 — Biz + SU psd fig 
onde e 
fy=D, (= Blhs)* = exp (—BUi;) — 1 (6.16) 
h=1 


Considerando que a energia entre as partículas é muito menor 


que a energia térmica, ou seja, Ui; €< kBT ou BUi;; < 1, conse- 
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quentemente, fi; < 1. Esta aproximação é justificável para um 
gás real em temperatura alta e/ou um gás muito diluído, pois a 
distância entre as moléculas é grande e, portanto, a energia poten- 
cial entre estas é pequena, como pode ser visto na fig. 6.1 . Como 
estamos querendo calcular o primeiro coeficiente do virial, é ra- 
zoável considerar baixa a concentração do gás (gás muito diluído) 
e, sendo assim, podemos fazer a seguinte aproximação 

[exp (— Bis) = Ia +fy)=1+ Bo dare (6.17) 

i<j i<j i<j 


De (6.17) em (6.15), obtemos 


E [711 / en) RES A 
l 


i<j 
= ve (x / en) De melo 

| i<j 
E 


h 


pois, o volume ocupado pelo gás é V = [ d3r;. Note que, 


Rn [7x1 / 4) Dt 
) i<J 
1 


- dm fem [ rulfat fia dec + noi) 


1 
= a fem ferra (fem fer feio) + 
A O (fam fem [er) Eos 
Tr, 7 ed Ed — 
o [diva | intros N (fam fem [erws) 
1 - 


1 = 4 o 
pm | Sina | divtao NNE 


1 
= SN (Ne D fm [dt 


8 
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Na última igualdade usamos o fato de que as partículas são idên- 


ticas e, portanto, 


fem fénto= [er fer 


para qualquer à < j. Por outro lado, o número de combinações 


possíveis de N partículas duas a duas é 


ou seja, 


fio,f13,..., ÍN-1 N 
2 — 


3 N(N-1) termos 


que corresponde ao número de integrais iguais somadas. 


assim, 


V—2 
Zi ES + N(N 1) fer [rifa 
Nm 
1 


Sendo 


No desenvolvimento da integral 1) usaremos o sistema de coorde- 


nadas ilustrado na fig. 6.2 . Logo, 


fem férato 


1 


= [em (ar 2 fia(r) [ao f acoso) 


I 
ca 
a, 

(94) 
aq 
3 


onde 
(DT) = tr [ dr r2 fio(r), 
0 


o qual não depende da posição da partícula 1. Portanto, 


Db a Va. 


(6.18) 
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Za 


a 
Dl ap 


W 
ma 
» 


- 
a SR 
— 


Xa YZ = 


Í 
Eseço 


Figura 6.2: Posição da partícula 2 em relação à partícula 1. (Figura 


adaptada da ref. [2].) 


Perceba que, em geral, fi; depende da temperatura [eq. (6.16)]. 
Sendo assim, 7 também é uma função da temperatura e, por isso, 


evidenciamos esta dependência. Com isso, podemos escrever 


N(N-1) 


Zn 2 1 
int + 2V 


n. 


Podemos considerar N(N — 1) = Nº, pois N > 1. Assim, temos 


Ainda podemos recorrer à seguinte aproximação 


1+ax)” =1+na, para x < 1. 
L 
080, Nm a 
Substituindo as egs. (6.14) e (6.19) na eg. (6.13), obtemos 
Ee V(CeamksT3/2 1 + ND) (6.20) 
NM! a 2V 


correspondendo à função de grande partição de um gás real com 
correções de primeira ordem em relação ao gás ideal. Perceba 


que se a interação entre as partículas for desprezada, fio = 0 = 


(116, 
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n = O resultando na função de grande partição do gás ideal, como 
esperado. 
Vamos agora, calcular a equação de estado para um gás real. 


Para isso, façamos 


NT 
nZ=z-nN!+NhnV+NIn(2amksT)2 + Nln [+ SH 


Considerando a aproximação de Stirling, nN! =x NnN —- Ne 
que o gás é muito diluído, ou seja, Ny(T)/V < 1, podemos usar 


a aproximação 
In(l+2)= a, parar <1. 


Com isso, temos 


nZ=NhV+ ia, a gre | 
Assim, 
av DP? = 7 B Nao (6.21) 
e 
SpZ= Nip tr or (6.22) 


Podemos calcular a pressão substituindo a eq. (6.21) na eq. 


(6.12). Com isso, temos 


pa NkBT NºkBT(T) 
na 2V?2 


(6.23) 


Comparando esta equação com a eq. (6.5), concluímos que o se- 


gundo coeficiente do virial é 
1 
A(T) = =otBTuD), (6.24) 


Perceba que a natureza da interação entre as partículas do gás 


está embutida em n(T). Portanto, além de A(T) depender da 
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temperatura, este também depende do tipo de interação entre as 
moléculas do gás. 

Analogamente, podemos obter a energia interna do gás substi- 
tuindo a eq. (6.21) na eq. (6.11). Logo, 


Nºkoto on 


3 
DR A ART 


(6.25) 


Exemplo 6.3.1. Calcule o segundo coeficiente do virial no caso 
de um gás composto de esferas rígidas não interagentes de raio ro 
e expresse a pressão deste gás com correção de primeira ordem. 


(Exemplo adaptado da ref. [2].) 


Solução: Para moléculas impenetráveis e não interagentes, temos 


00, se r<2ro; 
Uiol(r) = (6.26) 
0, se r>2ro, 


onde r é a distância entre os centros das moléculas. 


Sendo assim, temos 


Uio(r —l, se r<2ro; 
fialr) = exp | | 1= (6.27) 
B O, se r>2r o. 
Portanto, 
0,0) 
= dg dr r2 fy(r) 
0 
2ro (9,0) 
= 47 dr rº fr) + dr r2 fyo(r) 
ro Nm om” 
—- 0 
2 
= e p3 [É"0 = 327r% 
0 3 
Consequentemente, 
167rê 1 
MesbaT O|)=> kpTb 
É ( 3 ER 
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onde N4 é o número de Avogadro e 
1 3 
b=N A ( o) A 


corresponde ao quádruplo do volume de um mol de molécu- 


las. Com isso a equação de estado das variáveis de estado P, 


VeT em primeira ordem de aproximação é 


pa NkBT | NºkpTb 
oo TC veNa 


ou, como R = Nakp (constante universal dos gases), temos 


O NRT  NºRTb 


P 
V ve 


Perceba que nossa aproximação não é suficiente para obter 
o primeiro termo da correção da energia interna do gás em 
relação ao gás ideal, pois neste exemplo 7 não depende da 
temperatura [ver eq. (6.25)]. Para isso, seria necessário con- 


siderar mais termos na aproximação da eg. (6.17). 


6.4 Conclusão 


Nesta aula, generalizamos o formalismo de mecânica estatís- 
tica clássica abordado na aula passada para tratarmos sistemas 
partículas interagentes. Com isso, estudamos gases reais de forma 
aproximada. Vimos que expandindo a interação entre os constitu- 
intes do gás, somos capazes de encontrar correções para a aprox- 
imação das equações de estado do gás real em relação às do gás 


ideal. 
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6.5 Resumo 


Quando um gás real (partículas interagentes) possui uma con- 
centração (número de partículas por volume) baixa (gás muito 
diluído), podemos tratá-lo aproximadamente como um gás ideal. 
Por exemplo, a equação de estado das variáveis P, V,e T pode ser 
escrita como uma expansão polinomial em N/V (concentração de 


partículas) da seguinte forma: 


NRO a AN NC AR 
E To qo tg os 


P 
onde 4, Be € são chamados de coeficientes do virial. 
Para tratar estatisticamente um gás real, sofisticamos o formal- 


ismo de mecânica estatística visto na aula passada, onde apresen- 


tamos a função de grande partição 
Z = ZuZim 


onde Zcy é a função de grande partição do gás ideal, que é 
N 


V 
Zai — a (mmksT) NA, 


Zin = [7511 /n) [es (su), 


i<j 


Ui; =Ulri,T;). 


Com a função de grande partição podemos obter P e U da seguinte 
forma: 


Õ 
= 2 
U=kp7 2T In Z, 
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(6) 


Calculamos Zi;nt com aproximação de primeira ordem e obte- 


mos 


onde 
(TD) = tr [ dr r2 fio(r), 
0 


r=|r79-—m|éa distância entre as partículas 1e2,e 
fo(r) = exp [-BUho(r)] — 1. 


Com esta aproximação, encontramos 


pa NkBT NkBT(T) 
o A 9V2. 


A(T) = —5ksTn(T) 


NºkpTº On 


3 
NE Ned ES 
5 MEL ay BT 


6.6 Atividades 


ATIV. 6.1. Calcule a primeira correção em relação ao gás ideal 
para a pressão e para a energia interna de um gás composto de 
esferas impenetráveis e fracamente interagentes, onde há uma pe- 
quena repulsão entre estas quando não se tocam. (Atividade adap- 


tada da ref. [2].) 


Comentário: Este problema é análogo ao resolvido no exemplo 


6.3.1. No entanto, aqui as moléculas possuem uma interação 


(1205 
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muito pequena para r > 2r9, o que permite uma expansão 
em exp(— Bjo) até primeira ordem. Uma consulta na ref. 


[2] poderá auxiliar na resolução. 


ATIV. 6.2. Uma equação de estado para descrição de gases reais 


proposta empiricamente por van der Waals é 


(» + 2) (VN) = NRT. 


Escreva esta equação na forma do virial e compare com o resultado 


da atividade anterior. 


Comentário: A forma do virial equivale à expansão polinomial 


vista da eg. (6.5). 


ATIV. 6.3. Usando o resultado da ativ. 6.1 , calcule o trabalho 
realizado por um gás quando este se expande isotermicamente de 
Vi a VW. Aplique o resultado para um mol hidrogênio a 300 K 
quando este se expande de um volume de 3 x 10-2m? a um volume 
de 5 x 10-2m?. Compare com o valor obtido usando a expressão 


do gás ideal. (Atividade adaptada da ref. |2].) 


Comentário: Lembre-se do trabalho realizado por um gás visto 


na sec. 1.4.1 . Use os dados da tab. 6.1. 


ATIV. 6.4. A temperatura de Boyle de um gás real é a qual torna 


nulo o segundo coeficiente do virial. Mostre que esta temperatura é 
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Substância a (10-82) b (10525) 
Hélio 3,457 0,02370 
Hidrogênio 24,68 0,02661 
Neônio 21,35 0,01709 
Nitrogênio 140,8 0,03913 
Oxigênio 137,8 0,03183 
Amônia 422,5 0,03707 
Dióxido de Carbono 364,0 0,04267 
Água 553,6 0,03049 


Tabela 6.1: Constantes de van der Waals. 


iguala a/Rb. Calcule esta temperatura para todos os gases listados 


na tab. 6.1. (Atividade adaptada da ref. [2].) 


Comentário: Use o segundo coeficiente do virial para o gás de 


van der Waals. 


6.7 Próxima aula 


Na próxima aula estudaremos o calor específico de gases ideais 
mono e poliatômicos. Além disso, apresentaremos o princípio de 


equipartição de energia. 
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